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Abstract 

The thermodynamic properties of the quark gluon plasma (QGP) as well as its phase diagram 
are calculated as a function of baryon density (chemical potential) and temperature. The QGP 
is assumed to be composed of the light quarks only, i.e. the up and the down quarks, which 
interact weakly and the gluons which are treated as they are free. The interaction between quarks 
is considered in the framework of the one gluon exchange model which is obtained from the Fermi 
liquid picture. The bag model is used, with fixed bag pressure (B) for the nonperturbative part and 
the quantum chromodynamics (QCD) coupling is assumed to be constant i.e. no dependence on 
the temperature or the baryon density. The effect of weakly interacting quarks on the QGP phase 
diagram are shown and discussed. It is demonstrated that the one gluon exchange interaction for 
the massless quarks has considerable effect on the QGP phase diagram and it causes the system 
to reach to the confined phase at the smaller baryon densities and temperatures. The pressure of 
excluded volume hadron gas model is also used to find the transition phase-diagram. Our results 
depend on the values of bag pressure and the QCD coupling constant which the latter does not have 
a dramatic effect on our calculations . Finally, we compare our results with the thermodynamic 
properties of strange quark matter and the lattice QCD prediction for the QGP transition critical 
temperature. 
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I. INTRODUCTION 



The quark gluon plasma (QGP) usually is denned as the phase of Quantum Chromody- 
namics (QCD) in which the quarks and gluons degrees of freedom, that is normally confined 
within the hadrons, are mostly liberated. The possible phases of the QCD and the precise 
locations of critical boundaries or points are currently being actively studied. In fact, re- 
vealing the QCD phase transition structure is one of the central aims of the ongoing and 
the future theoretical and the experimental research in the field of the hot and/or the dense 
QCD [IHl]. It is about thirty years since the study of the hot and the dense nuclear matter 
in the form of the QGP has been started. The experiments at the CERN's super proton 
synchrotron (SPS) first tried to create the QGP in the 1980s and 1990s: The first hints of 
the formation of a new state of matter was obtained from the SPS data in terms of the 
global observable, the event-by-event fluctuations, the direct photons and the di-leptons. 
The current experiments at the Brookhaven national laboratory's relativistic heavy ion col- 
lider (RHIC) are still continuing these efforts. In the April 2005, the formation of the quark 
matter was tentatively confirmed by the results obtained at the RHIC . The consensus of 
the four RHIC research groups was in favor of the creation of the quark-gluon liquid at the 
very low viscosity [5HH]. 

Since, in this new phase of matter, the quarks and the gluons are in the asymptotic 
freedom region, one expects that they interact weakly. So, the perturbative methods can be 
used for such a system. The asymptotic freedom suggests two procedures for the creation 
of the QGP: (i) The recipe for the QGP at high temperature. If one treats the quarks and 
the gluons as the massless noninteracting gas of molecules, such that the baryon density 
vanishes, the critical temperature which above that, the hadronic system dissolves into a 
system of quarks and gluons (QGP) is T c = 140MeV |9J. However, the modern lattice 
QCD calculation estimates the critical temperature, %, to be about 170 MeV [TU] (ii) The 
recipe for the QGP at high baryon density. At zero temperature, the critical baryon density 
required the transition to take place in % ~ 0.7/m~ 3 [9] i.e. four times the empirical nuclear 
matter density. On these grounds, one should expect to find the QGP in two places in the 
nature: Firstly, in the early universe, about 10 _5 s after the cosmic big bang or secondly, at 
the core of super-dense stars such as the neutron and the quark stars. This new phase of 
matter can also be created in the initial stage of the little big bang by means of the relativistic 
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nucleus-nucleus collisions in the heavy ion accelerators [5HE]. 

The critical temperature (critical baryon density) at the zero baryon density (zero tem- 
perature) has been obtained simply for the non-interacting, massless up and down quarks 
and gluons in the reference [9]. Some primary works in the zero baryon density and in the 
framework of the bag model have been also presented in the reference [H] (and the reference 
therein). For the more complicated field theoretical approaches see the references [T2l [13] 
and the references therein. So it would be interesting to perform a similar calculation to 
reference (9], but with the non zero interaction, for the region in which both T ^ and 
fis / 0. So, in this work we generalize above calculations for the region with the finite 
temperature and the baryon density by considering the weakly interacting quarks in the 
framework of the one gluon exchange scheme. Since it is intended to use the perturbative 
method in our calculation, it is assumed that a c < 1. The dependence of a c on the temper- 
ature and the baryon density is ignored [HI [15]. We perform calculations beyond the zero 
hadronic pressure approximation of reference [9] and the pressure of the excluded volume 
hadron [TEHT5] model is also taken into account to find the corresponding QGP transition 
phase diagram. 

So, the paper is organized as follows. In the section // we review the derivation of the one 
gluon exchange interaction formulas [20] based on the Landau Fermi liquid picture [2TI 122] . 
The section III is devoted to the calculation of the thermodynamics properties of the QGP. 
The result and its comparison with the strange quark matter [13] are given in the section 
IV. Finally a conclusion and summary are presented in the section V. 



II. THE LANDAU FERMI LIQUID MODEL AND THE ONE GLUON EX- 
CHANGE INTERACTION 

The Landau Fermi liquid theory describes the relativistic systems such as: the nuclear 
matter under the extreme conditions, the quark matter, the quark gluon plasma, and other 
relativistic plasmas. The basic framework of the Landau theory of relativistic Fermi liquids 
is given by Baym and Chin [21] • In this framework, one can evaluate the energy density of 
a weakly interacting quarks in QGP by the following formulas: 

£ &kin £poti (1) 
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with 

e V n P ^27r) 3 ' ^ 
£• 1 / runpol 

g g d 3 p g g d 3 k 

tpot - 2 y J M n P n *( 2 7r) 3 (27r) 3 ' 1 J 

where g g = g spin x g co i or = 2 x 3 = 6, is the degeneracy of quarks, e^n, £ P ot an d nj are 
the kinetic and the potential energies and the familiar Fermi-Dirac liquid distributions of 
our quasi-particles, respectively. In the equation (3), f p n k p ° l is the Landau- Fermi interaction 
function, which is a criterion for the interaction between the two quarks (the "unpoF su- 
perscript refers to the unpolarized quark matter [TU [15]). It is related to the two-particle 
forward scattering amplitude i.e.: 

runpol m q m q KA ( A \ 

Jp,k ~ -Mp,fc- l 4 J 

A4 Pi k is the commonly defined Lorentz invariant matrix element. Now, using the Feynman 
rules for the QCD, Ai Pt k can be calculated. Since, the direct term is proportional to the 
trace of the Gelman matrices, the color symmetric matrix element is only given by exchange 
contribution: 

Mp, k = -G 2l -tr(^)x 

u{k)-/ u u{p)u{p)Yu{k) — _^ 2 , (5) 

where Q 2 = 47ia c (our choice is a c = 0.2). Since our system is unpolarized, it is possible to 
sum over all the spin states and get the average of this quantity as: 

2 g 2 1m\ - k.p 

Mp ' k =9^ q (k-p) 2 ■ (6) 
Then for massless up and down quarks, 

G 2 1 

Jp,k ~ "7T l„l I U • \') 



So, the total energy of interacting quarks can be evaluated by using the equations (1), (2), 
(3) and (7). In the next section, the various thermodynamic quantities of the interacting 
QGP like the internal energy, the pressure, the entropy, etc are calculated. 
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III. THE THERMODYNAMIC PROPERTIES OF THE QGP 



We begin with calculation of the free ultra-relativistic massless quarks partition function 
density, which follows by its non-interacting gluon contribution's. Since we are interested in 
the heavy ion processes, it is assumed that the initial state to be the nuclear matter with 
the equal neutron and proton densities So, 

n u = n d , (8) 

and the contributions of u and d quarks are equal. Then write the Fermionic partition 
function for each quarks is written as (V is the volume), 

= J [l n (l + e -^-»)) + l n (i + e -M+*))}^Ktf dk . ( 9 ) 
This integral can be evaluated analytically i.e., 

V 3(2tt)3 P 1 60 + 2 + 4 h 1 ' 

where v = nf3 and j3 — ^ (in the Boltzman factor unit). 

Now we are in the position to evaluate the various thermodynamic quantities, such as the 
energy density, the entropy and the free energy from the above partition function, especially 
we have, 

where G q = 2g q is the quark degeneracy for two, up and down quarks, flavors. All of 
the thermodynamic quantities are obtained as a function of chemical potential (/i) and 
temperature (T). The temperature is a suitable experimental quantity but the chemical 
potential is not. So, it is better to rewrite the thermodynamic quantities as a function of 
the temperature and the baryon density, instead of the chemical potential. Since the baryon 
number of a quark is |, we have, n = |n g , and: 

n=^(7r 2 T 2 /i + /i 3 ). (12) 

Then, the chemical potential can be found as a function of n and T, from the above equation 
i.e., 




/97T 4 „ 7T 6 T 6 V / 37T 2 /97T 4 . 7T 6 T 6 . 




FIG. 1: The density dependence of the free energy per baryon at two different temperatures. 
The solid and dash curves are for the weakly interacting QGP and the strange quark matter, 
respectively. 

So our thermodynamic quantities become a function of baryon density and temperature. 
For the energy density of gluons we have, 

r k A-ng g k 2 dk 
£ 9 = J e k/r _ I ( 2vr )3 

= ft£r>, (14) 

where again g g = 16 is the degeneracy of gluons. 

Our system is ultra-relativistic, so, there is a simple relation between the pressure and 
the energy density: 

V = \e. (15) 
The bag pressure and the vacuum energy contributions should be included, in addition to 
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FIG. 2: The temperature dependence of the free energy per baryon at two different baryon densities. 
The solid and dash curves are for the weakly interacting QGP and the strange quark matter, 
respectively. 



the quark and the gluon contribution energies, 

£QGp(n,T) = e q (n,T) + e g (T) + B 
V QGP (n,T) = V q (n,T)+V g (T)-B. 

Having the pressure, the entropy density of system is evaluated, 



(16) 
(17) 



(18) 



To study the phases of the QGP, one should concentrate on the QGP and the hadronic 
pressures, V QGP and V Ha dron [TSHU], respectively, i.e. for V QGP < V Ha dron, the system 
is in the confinement phase and the quarks and the gluons are inside the bag, but for 
"Pqcp > 'Pnadron, the quarks and the gluons pressures can overcome to the bag and the 
hadronic pressures and the system is in the de-confinement phase, i.e. the QGP phase is 
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FIG. 3: The density dependence of the entropy per baryon at two different temperatures. The solid 
and dash curves are for the weakly interacting QGP and the strange quark matter, respectively. 

created. So, the phase diagram can be extracted by solving the following equation: 

V QG p(n, T) = P q (n, T) + V g (T) - B — P H adr a n.(n(jM), T), (19) 

where in this work, both Puadrons = [9] and Puadrons ^ [I6HT8] cases are considered . 
For Puadrons 7^ 0, the excluded volume effect for the nuclear matter equation of state is used 
[T6HT8] to calculate the hadronic pressure (equation (29) of reference [17] ) i.e., 

'pideal ( . n~\ 

PHadrvnsHri, T) = P Had rons{», T) « f ^ons^ ' ) (2Q) 

[1 + vn{fi)\ 

where PH^ r<m8 {fJ',T) is the pressure of free ideal nucleonic (Fermion) matter [181 EE] (here, 
degeneracy is 4, m ~ 939MeV, v = 4|7rr 3 and r ~ 1.2 fm, e.g see the equation (9) of the 
reference j!8j). 

Now, the interaction energy between the quarks can be calculated by using the Landau 
Fermi liquid model and the results derived by us in the previous section. The potential energy 
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FIG. 4: The entropy per baryon as a function of temperature at two different baryon densities. 
The solid and dash curves are for the weakly interacting QGP and the strange quark matter, 
respectively. 

density of the interacting massless up or down quarks is found by using the equations (3) 
and (7), 



pot(^T) = - / -G -jn v n h , n _ Na , _^ , (21) 



1 f \ c2 \_ g q d 3 pg q d 3 k 

2 J 9 y pk npUk (2tt)3 (2tt) 3 

where g q = 6 is the quark degeneracy of one flavor. By using the Fermi-Dirac distribution 
and the value of each quark flavor potential energy i.e. the equation (21), an analytical 
formula for the total potential energy density of the two quark flavours is written as follows: 

'POt \ 



Note that, the potential energy for massless quarks is always positive. So, the interaction 
between quarks inside the bag is repulsive and it helps the interacting quarks and the gluons 
to penetrate from the bag more easily, rather than the noninteracting case, and further more, 
the one gluon exchange interaction, because of its repulsive properties, makes the conditions 
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FIG. 5: The QGP equation of state as a function of baryon density at two fixed temperatures. 
The solid and dash curves are for the weakly interacting QGP and the strange quark matter, 
respectively, but without the bag constant. 

easier for the system to make the transition to the QGP phase. 

The internal energy density of the QGP is evaluated by performing the summation over 
the interacting and the non-interacting parts of the energy density of quarks, the vacuum 
energy and the gluon energy density which were calculated before, and having that, the 
other thermodynamic quantities of the QGP are found. 

As it was pointed out before, it is assumed that the bag (hadronic) pressure is the one 
has been used in the reference [9] i.e. B = 208MeVfm~ 3 (zero), in order to compare our 
phase diagram with this reference, and since it is intended to compare our results with the 
reference [H], the QCD coupling constant os chosen to be a c = 0.2. On the other hand, for 
the none zero hadronic pressure, the baryonic chemical potential and the bag pressure are 
varied to reach to the critical temperature of the lattice QCD predictions [ID] , i.e. T c = 170 
MeV. 
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FIG. 6: The pressure at zero temperature for different coupling constants as a function of baryon 
density. 

IV. THE RESULTS OF THE QGP CALCULATION 

We begin by presenting the calculated free energy per baryon for both the QGP and the 
strange quark matter [13] in the figure 1 (2), as a function of baryon density (temperature) 
at the two different temperatures (baryon densities). The free energy for the QGP is larger 
than those of strange quark matter, since we know that the strange quark matter should be 
more stable than the QGP and therefore the strange quark matter free energy should be 
smaller than that of the QGP. As one should expect, the free energy increases (decreases) 
by increasing the baryon density (temperature). While, the QGP has less temperature 
dependent with respect to the strange quark matter, they have similar density dependent 
at fixed temperature. 

Similar comparisons are made for the entropies in the figures 3 and 4. The entropy per 
baryon for the QGP is an increasing (decreasing) function of temperature (baryon density) 
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FIG. 7: The phase diagram for the interacting and the noninteracting QGP. 

and it is smaller than that of the strange quark matter [13] . Again their dependence on the 
density is the same, but they behave especially differently at larger temperatures. 

The plots of the equation of states of both the QGP (without the effect of constant bag 
pressure) and the strange quark matter as a function of baryon density at the two different 
temperatures are given in the figure 5. The QGP equation of state is much harder than 
that of strange quark matter at the same baryon density and temperature. 

The pressure of weakly interacting QGP for the two different QCD coupling constants, 
and the noninteracting QGP (without the effects of constant bag pressure) at zero tempera- 
ture as a function of baryon density is shown in the figure 6. The increase in the interaction 
strength makes the pressure to rise, and therefore at the smaller baryon densities, the pres- 
sure of quarks becomes equal to the bag pressure. So, the interaction facilitates the quarks 
transition to the deconfined phase at lower density. The QCD coupling constant also plays 
the same role, i.e it will reduce the transition density. 

Finally, the phase diagrams for both the interacting and the noninteracting QGP are 
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FIG. 8: The phase diagram for the interacting QGP by taking into account the nuclear matter 
excluded volume pressure at four different bag pressures (MeV). 

shown in the figure 7 for VHadrons(n(fi),T) = 0. The one gluon exchange interaction which 
is repulsive, causes to get the QGP at the smaller baryon densities and temperatures. As it 
was pointed out before, the reason is very simple, the repulsive interaction between quarks 
helps them to escape from the bags. So the formation of the QGP happens much easier 
for the interacting quarks than the noninteracting one. But, the critical temperature is 
about lAOMev, which is much less than the lattice QCD suggestion of 170MeV. In the 



figure 8, the hadronic pressure has been also taken into the account (see the equations (19) 



and the reference 20). The slashed area is the forbidden region, i.e. for the bag pressure 



approximately larger than 200MeV there is no critical temperature with the zero baryonic 
chemical potential (density). With the bag pressure about 441MeV (note that the bag 
pressure estimated to as large as 500 MeV jl]) and the low baryonic chemical potential 
(density) it is possible to get results near to the lattice QCD prediction, i.e. % = 170 MeV. 
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V. CONCLUSION AND SUMMARY 



In conclusion, the one gluon exchange interaction was used to evaluate the strength of 
potential energy of the QGP in the Fermi liquid model. By calculating the QGP partition 
function, the different thermodynamic properties of the QGP as a function of baryon density 
and temperature for the both interacting and the non-interacting cases were discussed. It 
was found that the QGP internal and free energies are much larger than those of strange 
quark matter. On the other side, if we consider the massive quarks like the strange quarks 
in our the QGP, the potential energy becomes a negative quantity, but for the massless 
quarks it is always a positive quantity, therefore the internal and the free energy densities 
of strange quark matter become smaller than the QGP ones. We have seen how the one 
gluon exchange interaction for the massless quarks affects the phase diagram of the QGP 
and causes the system to reach to the deconfined phase at the smaller baryon densities and 
temperatures. Our results depend on the values of bag pressure and the QCD coupling 
constant which the latter does not have a dramatic effect on our results. The increase of the 
hadronic and the bag pressure can improve our results toward the lattice QCD calculations. 
In the future works, we could adjust our phase diagram to get the relation between the bag 
constant and the QCD coupling constant. On the other hand it is possible to generalize our 
method for the non-constant QCD coupling and the bag pressure. Finally, we can also add 
the interaction between the gluons to our present calculations. 
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